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INTRODUCTION: DMS and Cloud Formation in the Remote Marine Boundary Layer 


The general course of events in the formation of a marine cloud begins with the emission 
of species which can eventually serve as nuclei around which water can condense to form a cloud 
droplet. In remote marine regions, cloud condensation nuclei (CCN) are primarily composed of 
sulfate, in either its acid or ammonium salt form. Most sulfate in these regions is the product ot 
atmospheric oxidation of dimethyl sulfide (DMS), a reduced sulfur gas that is released by 
phytoplankton at the ocean surface. Therefore, in order to effectively quantify the links in the 
cloud-formation cycle, one must begin with a well-defined description of the atmospheric 
chemistry of DMS. The intent of this project has been to initiate development of a 
comprehensive model of the chemistry and dynamics responsible toi the formation ot clouds in 
the remote marine boundary layer. The primary tool in this work has been the Global/Regional 
Atmospheric Chemistry Event Simulator (GRACES), a global atmospheric chemistry model, 
which is under development within the Atmospheric Chemistry and Dynamics Branch ot NASA- 
Ames Research Center. In this effort, GRACES was used to explore the tiist chemical link 
between DMS and sulfate by modeling the diurnal variation of DMS. 

BACKGROUND: Previous Efforts to Model the DMS Diurnal Oxidation Cycle 

Examples of other efforts to model the diurnal variation of DMS are described in M.O. 
Andreae, et al. [1985], Koga and Tanaka [1993], Saltzman, et al. [1993], Yvon, et al. [submitted, 
1994], and Suhre and Rosset [1994], Recent efforts to model the rate of formation of CCN 
[Russell, et al. 1994; Suhre and Rosset, 1994], treat the initial oxidation of DMS as manifested 
by its diurnal cycle in a limited fashion, assuming a parameterization that does not account tor 
variations in flux and background atmospheric chemical composition. In this work, we consider 
these factors and the role they may play in dictating the rate of DMS oxidation. 

The following questions have been considered in this study: 1) How well can a zeio- 
dimensional chemical model predict the air concentrations of dimethylsulfide (DMSa) that were 
measured by the Andreae group [1994]? What are the implications, by extension, for higher 
dimensional DMS-CCN models which use zero-dimensional or simpler calculations as a first 
step in calculation? 2) Are there other significant sinks for DMS besides hydroxyl radical (OH) 9 
3) Is there evidence of a photoactive oxidant which becomes active at sunrise? 
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APPROACH: 


The Data 

The data considered in this study is composed of a 6-week continuous set ot 
measurements including hourly atmospheric and sea water concentrations of DMS. ozone, aii 
temperature, solar radiation, as well as halt-daily radiosonde measurements ot the maiine 
boundary layer height [Andreae, et al„ 1994]. DMS flux values were calculated using the 
measured wind speed, sea water temperature and the sea water DMS concentration. This data set 
was collected during a period of remarkably stable meteorological conditions, resulting in clear 
diurnal variation in hourly atmospheric DMS (DMSa) measurements. Analysis of the daily 
averages of the DMSa measurements versus the calculated DMS flux for February 12 through 
March 13 shows a strong correlation (r^ = 0.91). An average of the hourly atmospheric DMS 
measurements taken during the March 8 - 16 period have been used in a one-dimensional 
modeling study of sulfur chemistry and aerosol formation in the marine boundary layer [Suhre. et 

al., 1994]. 

The GRACES -DMS Diurnal Variation (GRACES-DDV) Model 

The primary oxidant of DMS is expected to be hydroxyl radical (OH). Nitrate radical 
(N0 3 ) is also thought to play a role in DMS removal at night. The concentrations of these and all 
other species are calculated by the GRACES model by integration of a full set of atmospheric 
reactions of importance in the troposphere. OH concentration is dependent upon the availability 
of ultraviolet light of the correct frequency to photolyze ozone. The photolysis rate of ozone and 
other photolyzible species is calculated using a 5-Eddington technique, modified for the presence 
of clouds, for an assigned surface albedo and for the solar radiation appropriate for the time ot 
year and global location being modeled. These photochemical reaction rates are calculated based 
on tabulated absorption and quantum yield data. In this study, the ocean surface albedo was set 
to 0.1. The rates for non-photochemical reactions are derived from DeMore, et al., JPL 
evaluation, Atkinson, et al. IUPAC evaluation, the Paulson and Seinfeld condensed isoprene 
mechanism, Lurmann et al. condensed mechanism for C 2 C 3 , C 6 -C 8 (aromatic) species. 
Modifications to the hydrocarbon mechanism to include hydroperoxide chemistry were made as 
described in Chatfield and Delany. Rates for oxidation of DMS by OH and N0 3 radicals are 
those reported by Hynes, et al. [1986] and Atkinson, et al. [1992], The full set of reactions is 
integrated by a differential-algebraic solver (DASSL) in fifteen minute time steps. 
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Chemical species are supplied to the model in two ways. All species are initially 
assigned a specific mixing ratio. Table 1. provides a list of the species and concentrations, as 
well as the temperature, standard altitude, gas number density and surface albedo used in a 
typical model calculation. 

Table 1. 


General Variables 

3.3E+3 cm 

Altitude 

Temperature 

298.1 K 

Number Density 

2.4E+19 (#/cm 3 ) 

Surface Albedo 

0.1 

Chemical Species 

Mixine Ratio 

H 2 0 

2.39E-02 

o 3 

1.82E-8 

NO 

0 

no 2 

5.00E-1 1 

no 3 

0 

hno 2 

0 

CO 

5.00E-8 

ch 4 

1.70E-6 

h 2 o 2 

7.00E-10 

HC1 

1.00E-10 

h 2 

5.02E-7 

DMS 

1.23E-10 

c 3 h 6 

3.60E-1 1 

c 2 h 4 

5.60E-1 1 

c 2 h 6 

5.60E-10 

c 3 h 8 

3.60E-10 

c 4 h 8 

5.00E-1 1 


The values for CH 4 , H 2 , and H 2 0 2 used are global averages. The mixing ratio for water was 
calculated from the daily dew points measured aboard the Meteor. Ozone (O 3 ) is maintained at 
its initial concentration by the model code to simulate its replenishment due to probable marine 
boundary dynamics. This assumption is made to fit observed concentrations. The ozone mixing 
ratio used was the average of the hourly measurements for the day being modeled. NOx and CO 
concentrations were estimated for the clean conditions observed during the cruise. Source terms 
for DMS and other reactive non-methane hydrocarbons are modeled as "appearance" rates, 
calculated from measured fluxes in the case of DMS or from measured background values for 
ethane (C 2 H 6 ), ethene (C 2 H 4 ), propane (C 3 H 8 ), propene (C 3 H 6 ) and butene (C 4 H 8 ). Here it is 
assumed that the background concentrations of other biogenic hydrocarbons are proportional to 
the DMS background value. Rudolph and Johnen, in 1990, measured light hydrocarbons in the 
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centra] Atlantic region crossed at the mid-point of the Meteor cruise. Andreae, et al. measured 
DMSa concentrations that were a factor of two smaller than those measured in the part of the 
cruise of primary interest in this study. Hence, a simple estimate was made that the background 
levels of the light non-methane hydrocarbons in the region considered here are double those 
measured by Rudolph and Johnen. Appearance rates were then calculated to maintain a steady 
state against removal by OH at its daily average concentration. Table 2. lists the non-methane 
hydrocarbon source terms included in the model and the base values from which they were 
estimated. 


Table 2. 

Chemical Species 

Ethene (C 2 H 4 ) 

Ethane (C 2 H 6 ) 

Propene (C 3 H 6 ) 

Propane (TNHg) 

Butene (C 4 H 8 ) 

1. Rudolph and Johnen. JGR 95(D12). p. 20583. November 20, 1990. 

2. Donahue and Prinn. 1990. 

To calculate a DMS appearance rate from the hourly flux measurements taken aboard the 
Meteor, the average of the radiosonde measurements which gave a marine boundary layer height 
of 700 m was used. Here, the assumption was made that the marine boundary level is 
instantaneously well-mixed. 

Since the principal focus of this work is the chemistry of DMS oxidation, the model is 
operated in its 0 -dimensional mode. It is assumed that the diurnal change in DMS concentration 
can be calculated from the difference between the source term and its loss by chemistry. In this 
study, DMS dry and wet depositional processes are considered unimportant, due to its extremely 
low solubility in sea and cloud water. Other loss processes such as vertical entrainment by 
clouds at the top of the marine boundary layer will be considered later. Depositional losses of 
other species are also not considered within the model at present. 


Low Productivity, 
Measured (ppt) 

25 1 
278 1 
18 1 

25 2 


Mod. Productivity, 
Estimated (ppt) 

50 

556 

36 

360 

50 


Mod. Source Rate, 

Calc. (#/cm3*s) 

6.5E+4 

2.5E+6 

1.6E+5 

5.9E+6 


Analysis of the Meteor Dataset 

As indicated above, the primary assumption in GRACES-DDV 0-dimensional chemical 
model is that the sole source of DMS is due to its flux from the ocean surface and its removal is 
due exclusively to oxidation. Thus, a correlation necessarily exists in the model between the rate 
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of flux and the concentration of DMS. In order to determine the impact of chemistry on the 
atmospheric DMS concentration, a data set which shows a strong correlation between flux and 
DMSa is required. Given the strong correlation between average daily DMS concentration and 
average daily DMS flux rate shown by Andreae, et al.. an effort was made to determine the 
degree to which the hourly flux values correlated with the hourly DMSa measurements foi 
individual days within their data set. The measured flux and DMSa values were plotted and a 
linear regression calculated for each day. The individual days which gave the best linear 
correlations, i.e. R > 0.8. were modeled, along with the subset chosen tor modeling by Suhre. et 
al [19941. 

Individual days were also normalized against their average DMS concentration, as a 
means of comparing the relative change in DMSa observed in the data, versus the change 
calculated in the model. 

Furthermore, because the chemistry of interest in this study is thought to occur at specific 
times within the solar day, a similar analysis of the flux versus DMSa for night, eaily moining. 
mid- and late day was carried out. Individual days were divided into these time regimes 
according to flux of ultraviolet radiation measured on the particular day. Night time was defined 
as the periods beginning at midnight until sunrise and one hour after sunset through the following 
midnight. Early morning was chosen as the first 2-3 hours starting at sunrise and 
ending at the midpoint of the rising edge of the UV flux curve. Midday was contained between 
the midpoints of the UV flux rise and fall. Late day was the interval between the end of the 
midday period and sunset. Table 3. shows the hours selected for each interval by date. 


Table 3. 
Date 

February 20 
February 28 
March 4 
March 8 
March 9 
March 10 
March 1 1 
March 12 
March 13 
March 14 
March 15 
March 16 


Night-time 

Earlv Morning 

0-7,21-23 

8-10 

0-5, 19-23 

5-8 

0-6, 19-23 

7-9 

0-5, 18-23 

6-7 

0-5, 18-23 

6-7 

0-5. 19-23 

6-7 

0-4, 18-23 

5-6 

0-5, 19-23 

6-7 

0-5, 19-23 

6-7 

0-5, 19-23 

6-7 

0-4, 18-23 

5-6 

0-4. 18-23 

5-6 


Middav 

Late Dav 

11-17 

18-20 

9-15 

16-18 

10-15 

16-18 

8-15 

16-17 

8-15 

16-17 

8-16 

17-18 

7-15 

16-17 

8-16 

17-18 

8-16 

17-18 

8-16 

17-18 

7-15 

16-17 

7-15 

16-17 
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The linear rate of change in DMS concentration at early morning and midday was also calculated 
explicitly for both the model and data for individual days. 


Model Calculations Undertaken 

A subset of daily measurements was modeled, using the 24 hour average flux, ozone, 
DMS level, water mixing ratio and air temperature for each day. The model was also provided 
with the latitude, longitude and date which the measurements were taken to assure that 
appropriate solar radiation flux values were used in the photochemical rate calculations. Flux to 
DMSa correlation plots were constructed and the change in early/midday DMS w-as calculated 
for comparison to the data. Table 4. shows the measured variables used in the model 
calculations. 


Table 4. 
Pate 


Feb 20 
Feb 28 
Mar 4 
Mar 8 
Mar 9 
Mar 10 
Mar 1 1 
Mar 12 
Mar 13 
Mar 14 
Mar 15 
Mar 16 


DMS 

DMSg 2 

Ozone 3 

Water 4 

Air 

Latitude 6 

Long. 6 

Flux 1 




Temp 5 



2.81 

32.35 

14.12 

2.30 

25.75 

-18.63 

-24.94 

3.50 

56.42 

13.35 

2.40 

25.00 

-22.82 

-14.91 

9.15 

155.6 

9.458 

2.38 

24.32 

-19.60 

-6.761 

8.90 

8.32 

138.5 

133.9 

16.33 

18.61 

2.32 

2.39 

24.12 

24.95 

-19.00 

-19.00 

-0.332 

1.34 

9.30 

146.9 

18.17 

2.27 

23.22 

-19.00 

3.00 

13.9 

193.8 

17.05 

2.26 

23.05 

-19.00 

4.72 

19.7 

167.3 

17.73 

2.18 

22.62 



24.3 

129.2 

18.47 

1.93 

22.03 



14.2 

99.22 

18.58 

1.77 

21.82 



13.6 

111.9 

19.17 

1.88 

21.95 



13.8 

139.2 

21.37 

1.82 

21.63 




Unit* 

1 . pmol/m 2 *day — mean measured value for given day 

2. parts per trillion — mean measured value for given day 

3. parts per billion — mean measured value for given day 

4. percentage by number — mean measured value for given day 

5. Celsius — mean measured value for the given day. 

6. degrees — location of the ship at noon on the given day 


General calculations were performed to assess the relative importance of each of the key 
variables in determining the calculated concentration of DMS. The sensitivity of OH 
concentration to CO concentrations between 40 and 60 ppb and NOx concentrations between 5 
and 50 ppt was evaluated. For a direct assessment of the consequences to the DMS 
concentration, the model was supplied with a +/- 20% range around the measured or estimated 
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values of DMS flux, ozone, water, carbon monoxide and nitrogen oxides. This calculation was 
especially useful in estimating the impact due to error in the guessed values of NOx and carbon 
monoxide, which were not measured aboard the Meteor. 

Given the interest shown in the atmospheric chemistry literature in the possibility ot a 
halogen-based, photochemically active oxidant, an evaluation of the behavior ot such an oxidant 
on the diurnal variation of DMS was carried out. A source term and the necessary photolysis rate 
information for hypochlorous acid (HOC1), a gas-phase species that might evolve from acidified 
sea salt aerosol, was included in the model. Due to the photochemical similaiity expected 
between HOC1 and HOBr, this evaluation is also relevant to HOBr. 



RESULTS 


Table 5. Correlation statistics for hourly DMSa vs. DMS flux values measured by day. 


Pate 

Julian Date 

February 12 

31819 

February 13 

31820 

February 14 

31821 

February 15 

31822 

February 16 

31823 

February 17 

31824 

February 18 

31825 

February 19 

31826 

February 20 * 

31827 

February 2 1 

31828 

February 22 

31829 

February 23 

31830 

February 24 

31831 

February 25 

31832 

February 26 

31833 

February 27 

31834 

February 28 * 

31835 

March 1 

31836 

March 2 

31837 

March 3 

31838 

March 4 * 

31839 

March 5 

31840 

March 6 

31841 

March 7 

31842 

March 8 * 

31843 

March 9 * 

31844 

March 1 0 * 

31845 

March 1 1 * 

31846 

March 1 2 * 

31847 

March 13 * 

31848 

March 14 * 

31849 

March 15 * 

31850 

March 16 * 

31851 

March 17 

31852 

March 18 

31853 

March 19 

31854 

March 20 

31855 

March 2 1 

31856 


Slope 

R Value 

-2.77 

0.436 

-7.22 

0.393 

6.26 

0.630 

13.7 

0.533 

10.3 

0.620 

4.23 

0.440 

0.602 

0.0232 

2.63 

0.337 

15.1 

0.895 

3.49 

0.129 

22.9 

0.740 

15.6 

0.721 

0.665 

0.0574 

-2.13 

0.126 

6.39 

0.330 

11.6 

0.384 

22.0 

0.816 

-1.67 

0.0345 

11.6 

0.688 

8.50 

0.477 

10.8 

0.831 

-2.80 

0.281 

17.41 

0.794 

5.56 

0.322 

8.56 

0.392 

8.96 

0.900 

-2.26 

0.159 

2.45 

0.336 

2.91 

0.248 

0.378 

0.129 

-0.852 

0.149 

0.679 

0.167 

-0.603 

0.0707 


In addition to the subset modeled by Suhre, et al. (March 8-16), February 20, February 28 
and March 4 were also chosen for calculation as having shown the best linear correlation 
between flux and DMSa. 
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Figures 1-12 show the results of the model calculations in relation to the hourly DMSa 
measurements. For comparison, DMS and UV flux are also shown. Figure 13 is a plot of the 
difference in DMS concentration, normalized for average background DMS, between the data 
and the model for the four most highly correlated days. Figure 14 shows the linear correlation 
between flux and DMSa for the data and the model for the same set of days. Figure 15 compares 
the model with the March 8-16 data set. Table 6. summarizes the amplitude, length and mean 
DMS values for the diurnal cycles modeled 


Table 6. Variablity in DMSa by day: Measured and Modeled 


Date 


Feb 20 
Feb 28 
March 4 
March 8 
March 9 
March 10 
March 1 1 
March 12 
March 13 
March 14 
March 15 
March 16 


Cycle Amp. 
Data 

1.1005 

1.0886 

0.77410 

0.67500 

0.61 100 

0.81850 

0.80090 

0.99800 

0.56790 

0.80500 

0.58900 

0.72590 


Cycle Amp. 
Model 

1.2200 

1.2572 

1.0500 

1.4370 

1.5257 

1.4938 

1.4357 

1.3228 

1.1792 

1.2542 

1.3287 

1.3874 


Cycle Leng. 
Data 

13.300 

1 1 .400 
12.770 
10.030 
9.0100 
9.9000 

11.350 
13.320 
10.310 
10.270 
8.8500 
9.5800 


Cycle Leng. 
Model 

6.8000 

7.2500 

7.4000 

7.1300 

6.8300 

6.7300 

7.1300 

6.4300 

6.5300 

6.9000 

6.7000 

6.7300 


Mean DMS 
Data 

32.350 

56.420 

155.59 

138.50 

133.94 

146.90 

193.80 

167.26 

129.20 
99.220 

111.89 

139.21 


Mean DMS 
Model 

32.940 

51.350 

160.80 

127.80 

110.36 
127.44 
203.67 
218.97 

262.36 
184.41 
175.93 
172.34 


Figure 16. compares the ratios of maximum measured to minimum measured DMSa to 
the same ratio calculated for the models of oxidation of DMS by OH alone, OH and N03 with 50 
ppt NOx and 100 ppt NOx present, respectively. The impact on the model trace by variation ot 
certain key variables is shown in Figures 17-19. Figure 20 shows the values for OH calculated 
by GRACES for a range of initial carbon .monoxide and nitrogen dioxide concentratons. Figure 
21 shows the partitioning of NOx within a 24-hour period within the model. Table 7 gives a 
comparison between the rates of change with time in DMS at midday (10AM-2PM) when the 
OH concentration is expected to be highest. These rates were calculated as linear fits to the data 
and model in this time interval. 
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R Value (Data) Slope (Model) R Value (Model) 


Table 7. 

Date Slope (Data) 


Feb 20 

-4.63 

Feb 28 

- 7.52 

March 4 

-4.68 

March 8 

-8.18 

March 9 

-14.1 

March 10 

-22.3 

March 1 1 

-13.0 

March 12 

-14.5 

March 13 

-7.49 

March 14 

-10.2 

March 15 

- 7.46 

March 16 

-22.1 


0.994 

-5.44 

0.97 1 

-11.7 

0.925 

-32.5 

0.684 

-31.9 

0.915 

- 27.5 

0.991 

-32.1 

0.986 

-52.1 

0.952 

-53.2 

0.775 

-61.1 

0.985 

-43.8 

0.921 

-43.0 

0.951 

-42.9 


0.967 

0.981 

0.991 

0.973 

0.963 

0.967 

0.975 

0.975 

0.981 

0.981 

0.976 

0.972 


Figures 22-25 show the results of the comparison of data to model for photochemically 
interesting time regimes. Parts a. are correlation plots for March 8-16. Parts b. are correlation 
plots for Feb 20, Feb 28, March 4 and March 9. 


Figures 26 and 27 show the results of a modeling of HOC1/HOC1 source terms that are 
25%, 50% and 100% of the DMS source strength. 


CONCLUSIONS, to date: 

The model consistently overpredicts the morning, and underpredicts the late afternoon 
concentration of DMS in all cases considered. The afternoon underpredication may be remedied 
by inclusion of a more complete sulfur oxidation mechanism than currently exists in the model. 
The consequence will be a larger sink for OH as the DMS oxidation products react further. At 
this juncture, it appears that an early morning drop off in measured DMSa that is not predicted by 
this chemical model. Due to the difficulty in finding a large set of days with a strong DMSa to 
flux correlation, it might be wise to consider the possibility of meteorological effects. The 
alternative suggestion is that a photochemically active oxidant accumulates at night, and then 
becomes active at sunrise. HOCl/HOBr behave photochemically in this manner, however, the 
source terms required to make a discernible difference in the modeled trace of DMS are so large 
as to be highly unrealistic. 
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Additional analysis is required to evaluate the days best modeled to determine the 
conditions under which a zero-dimensional model can be used most effectively to provide a 
DMS source for CCN. The time scales of chemistry and dynamics in the marine boundary layer 
must be factored into this analysis. It may be possible to adapt the GRACES model to allow foi 
hourly changes in DMS flux, which may impact the ability of the model to fit the early morning 
DMSa concentration, in particular. 
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Figure 13. 
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Figure 15a. 

Flux vs. DMSa: Measured and Modeled 
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Figure 15b. 

Flux vs. DMSa: Data and Model 
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Figure 15c. 

Flux vs. DMSa; Data and Model 
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Max/Min (DMSa Data) 


Figure 16. 

Ratio of Maximum [DMSa] to 
Minimum [DMSa]): Data and Model 


Max/Min (DMSa Data) 

■ • - - Max/Min (OH Only Calc) 

Max/Min (OH + N03, 50ppt NOx) 
a- - Max/Min (OH + N03, lOOppt NOx) 



Date 


29 




Trillion 


Figure 17. 








Figure 19. 
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Figure 21. 
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DMSa (ppt) 


Figure 22b. 


Night DMS Flux vs. DMSa: Data and Model 
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DMSa (ppt) 


Figure 23a. 

Measured vs. Modeled Early Morning DMSa vs. Flux 
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Figure 23b. 

Early Morning DMSa vs. Flux: Measured and Modeled 
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Figure 24b. 
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Midrday DMSa vs. Flux: Measured and Modeled 
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Figure 25a. 

Measured and Modeled Late Afternoon DMSa vs. Flux 
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Figure 25b. 


Late Day DMSa vs. Flux: Measured and Modeled 
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